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ABSTRACT 

Diffusive nuclear burning of H by an underlying material capable of capturing protons can readily 
consume H from the surface of neutron stars (NSs) during their early cooling history. In the absence 
of subsequent accretion, it will be depleted from the photosphere. We now extend diffusive nuclear 
burning to He, motivated by the recent observation by Ho & Heinke of a carbon atmosphere on 
the NS in the Cassiopeia A supernova remnant. We calculate the equilibrium structure of He on an 
underlying a capturing material, accounting for thermal, mass defect, and Coulomb corrections on the 
stratification of material with the same zeroth order \i e = A/Z. We show that Coulomb corrections 
dominate over thermal and mass defect corrections in the highly degenerate part of the envelope. We 
also show that the bulk of the He sits deep in the envelope rather than near the surface. Thus, even 
if the photospheric He abundance is low, the total He column could be substantially larger than the 
photospheric column, which may have implications for rapid surface evolution (w 1 yr timescales) of 
neutron stars. When nuclear reactions are taken into account, we hnd that for base temperatures 
> 1.6 x 10 8 K, He is readily captured onto C. As these high temperatures are present during the 
early stages of NS evolution, we expect that the primordial He is completely depleted from the NS 
surface like the case for primordial H. We also find that magnetic fields < 10 12 G do not affect our 
conclusions. Armed with the results of this work and our prior efforts, we expect that primordial 
H and He are depleted, and so any observed H or He on the surfaces of these NS must be due to 
subsequent accretion (with or without spallation). If this subsequent accretion can be prevented, the 
underlying mid-Z material would be exposed. 

Subject headings: diffusion - nuclear reactions, nucleosynthesis, abundances - pulsars: general - stars: 
abundances, interiors - stars: neutron 



1. INTRODUCTION 

The surface composition of neutron stars (NSs) re- 
mains an outstanding problem. Observations of young 
NSs have failed to find indisputable signatures of any 
particular element on their surface. Some NSs com- 
pletely lack any spectral features such as RX J1856.5- 
3754 (Burwitz et al. 2001, 2003; Pons et al. 2002; Ho et 
al 2007). However, when spectral features are seen, the 
degeneracies in our (currently) limited understanding of 
the chemistry and physics of highly magnetized material 
makes secure identification a challenge. Theory also fails 
to constrain the initial conditions, i.e., supernova fall- 
back, spallation, subsequent accretion, that would make 
a convincing case for a particular composition. 

Even after the initial composition is set, subsequent 
evolution can occur. For instance, the small mass (~ 
10~ 17 A/ Q ) and rapid diffusion time (w 1 s) at the NS 
photosphere suggest that it should consist of H due to 
gravitational settling. However, we showed in Chang 
& Bildsten (2003,2004; hereafter paper I and II) and 
Chang, Arras, & Bildsten (2004; hereafter paper III) that 
H is easily destroyed by diffusive nuclear burning (DNB). 
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The central idea behind DNB, which was first proposed 
by Chiu & Salpeter (1968; see also Rosen 1970), is that 
H can diffuse to great depth where the temperature and 
density are sufficiently large to consume H by proton 
captures onto heavier elements. In paper I and II, we 
showed that this process is so effective that we expect 
NS surfaces to be depleted of any primordial H. The ob- 
servation of H on the surfaces of NS would then point 
to late-time or continuous accretion. This conclusion is 
insensitive to the strength of the magnetic field and the 
size of an inert He buffer that sits between the H and the 
underlying proton capturing elements (paper III). 

With all the H consumed, one would expect surfaces of 
He. However, recent observations by Ho & Heinke (2009) 
suggest that the NS at the center of the Cassiopeia A su- 
pernova remnant has a carbon surface with an effective 
temperature T c = 1.8 x 10 6 K and radius R = 12 — 14 km. 
This paper shows that He is also vulnerable to diffusive 
nuclear burning on the surface of Cas A. All primordial 
H/Hc is consumed during its early cooling history, ex- 
posing the underlying material. If subsequent accretion 
does not cover this underlying material, we would expect 
a population of NSs with mid-Z surfaces. 

We first review the physics of DNB in $2j outlining 
the major results of papers I-III. We then discuss the 
physics of material in diffusive equilibrium on NSs in 
We show how to calculate the electric field in de- 
generate material accounting for thermal, mass defect, 
and Coulomb corrections. Using this electric field, we 
then calculate the structure of an envelope composed of 
He and C. We show that the diffusive tail of He pene- 
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trates sufficiently deep to be readily captured by C, and 
find that He is easily consumed in young, hot NSs for a 
range of a-capturing elements. In [J5l we conclude with a 
discussion of the implications of our work for young NS 
surface compositions. 

2. BASIC PHYSICS FOR HYDROGEN DEPLETION 

The diffusion rates on NSs imply that H will rapidly 
float to the photosphere where the density is p w 
1 g cm -3 and the temperature of T k 10 6 K. Hydrogen 
at this density and temperature will not burn in a Hub- 
ble time. However, 10 3 — 10 4 cm below the photosphere, 
the higher densities (p > 10 6 gcm _3 ) and temperatures 
(T ~ 1 - 10 X 10 8 K) allow H to burn. The H that diffuses 
down to this depth will be readily captured, driving a H 
current from the surface to the burning layer. Over time, 
this results in the total depletion of H. 

In paper I, we studied the case where the rate of H 
burning is sufficiently slow that H always remain in dif- 
fusive equilibrium, and the rate limiting step is the rate 
at which H is consumed by nuclear reactions in the burn- 
ing layer. This burning layer is defined by the competing 
effects of an exponential decline in H number density and 
a rising temperature. For cooler NSs (with effective tem- 
peratures T e < 10 6 K), DNB takes place in this nuclear- 
limited regime. For hotter NSs (T e > 1.5 x 10 6 K), DNB 
is not limited by nuclear reactions in the burning layer, 
but rather by the rate at which H can diffuse down to 
the burning layer, i.e., the diffusion limited regime (pa- 
per II). In this case, the burning occurs in a thin layer 
at which the nuclear burning time becomes comparable 
to the downward diffusion time, breaking the assump- 
tion of diffusive equilibrium (see papers I and II). Below 
this burning layer, the hydrogen concentration profile is 
rapidly cut off. 

In papers I and II, we demonstated that DNB can eas- 
ily consume all the H on NS surfaces early in their cooling 
history. We also studied the effect of magnetic fields on 
the rate of DNB and found that for fields < 10 12 G, 
this basic evolution is unaffected. In paper III, we ex- 
tended our study of DNB to magnetar surfaces. There, 
we showed that the high temperatures associated with 
magnetars, i.e., soft gamma repeaters and anomalous X- 
ray pulsars, allow H to be captured onto elements as 
heavy as Fe. In addition, we showed that the effect of an 
inert He buffer between the H and C does little to damp 
the effectiveness of DNB. Our calculation of the struc- 
ture of an (assumed) inert He on C in paper III was an 
early attempt at calculating the structure of a degener- 
ate He/C atmosphere (including Coulomb physics). We 
refine our methodology as we describe below. 

3. ELECTRIC FIELD IN DEGENERATE PLASMA 

We first discuss the compositional profile of degenerate 
material in diffusive equilibrium, set by the competition 
between gravity and the electric field. To simplify the 
discussion initially, we consider an ideal gas of electrons 
and ions. In section 3.2, we include non-ideal effects. 
The equations of hydrostatic balance of a plasma of ions 
and electrons in a plane parallel atmosphere arc 



dPi 

= -m (AiiTipg - ZieE) , 

dP e 

-7T— = -n e eE, 
or 



(1) 

(2) 



where ni (n e ) is the number density of ions (electrons), 
P. L (P e ) is the ion (electron) pressure, Ai and Zi are 
the mass and charge of the ion, g = GAf N s/i?NS is the 
gravitational acceleration, M-^s and -Rns are the mass 
and radius of the NS, and E is the electric field. Note 
that we ignore the electron mass, and we set the nuclear 
mass to Arrip, where m p is set to be the atomic mass 
unit. This small change, which differs from our defini- 
tions in papers I, II, and III, allows us to include mass 
defect corrections. However, by defining the atomic mass 
unit as m p rather than the more standard u, we main- 
tain consistent nomenclature with papers I, II, and III. 
The electric field in equation ([TJ and ([2]) results from the 
condition of charge neutrality n e — Ziiii. For an ideal, 
nondegenerate plasma (P e = ZPi), this electric field is 
eE = Aiiripg/ (Zi + 1). For a strongly degenerate plasma 
(P e ^> Pi), the electric field becomes 



Ai 

eE = —m p g. 



(3) 



If we now consider the forces on a trace ion species j 
immersed in a background ion species i, equation ([1} be- 
comes 



dr 



-Tijiripg I Aj — Zj 



A, 



(4) 



For Aj/Zj = Ai/Zi, the differentiating force between 
trace and background ion species is zero for the zeroth 
order electric field of equation ([3]). Note this is not 
the case of H on C in paper I, II, and III, which have 
very different p e = A/Z and so the zeroth order electric 
field imposes a differentiating force. Here ions with the 
same Ai/Zi experience no differentiating force (ignoring 
mass defect corrections) in a strongly degenerate plasma. 
Hence, higher order corrections to the electric field due to 
thermal effects and Coulomb interactions determine the 
differentiating forces between ion species and the equi- 
librium compositional structure. 

3.1. Thermal, Mass Defect, and Coulomb Corrections 

We now give a qualitative discussion of the thermal 
and Coulomb corrections to the electric field in degener- 
ate plasmas. We first consider corrections to the zeroth 
order electric field due to purely thermal effects because 
this physics is simpler and illustrative. Using equations 
HI) and ^ and the condition of charge neutrality, the 
electric field is (paper I), 



^ niZi (Aim p g + k B (dT/dz)) 

i 

{dP e /dT) (dPe/dUe)- 1 (dT/dz) 



EniZf 



+ n e (n e k B T/P e ) (dlnP e /dlnn e y 



-(5) 



where z = r — r re f is the height measured from a reference 
point, r rc f. Note that we have ignored the mass of the 
electron in the derivation of equation ([5] (see paper I for 
a detailed derivation). As expected, we reproduce the 
zeroth order electric field (eq. [3J ) from equation (JSJ) for a 
degenerate plasma. 

Now consider a two component plasma where n\ n-i 
(1 is the background and 2 is the trace) with the same 
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(X e , i.e., A\jZ\ = A<ilZ<i- For a degenerate plasma 
(P e ^ n e kBT), we expand equation ([5]) in the ra- 
tio of the thermal pressure to the degenerate pressure, 
n e ,ikBT/P e , to find 



eE = -m p9 \1 - _ 



(6) 



where 7 = 91nP e /91nn e , demonstrating that the ther- 
mal correction yields a differentiating force between ions 
with the same Ai/Zi. We should note that Hameury, 
Heyvaerts and Bonazzola (1983) derived a similar results. 

Plugging equation ([6]) into hydrostatic balance for each 
ion, we find 



9 In Pi 
dr 
<91nP 2 _ 
dr 



P9 

Zi pg 



(7) 



where we presume that the ions obey the ideal gas equa- 
tion of state, Pi = niksT. Subtracting these equa- 
tions from each other and using hydrostatic balance 
dP/dr = —pg, where P ~ P e , we find a simple power 
law between the concentration, nijn\, and pressure, P, 



d In nijn\ 
dhiP 



7 



£2 
Zi 



1 



(8) 



Equation © demonstrates that thermal corrections to 
the electric field yield a differentiating force which sepa- 
rates lighter ions from heavier ions with the same Ai/Zi 
in a degenerate plasma. 

Our discussion of purely thermal effects above, while 
illustrative ignores mass defect corrections and Coulomb 
corrections. Mass defect corrections arise from that fact 
that Aj/Zj is only « Ai/Zi for different i and j due 
the the slight variation in atomic mass. While this is 
usually ignored, as we have done in paper I, II, and III, 
the anonymous referee has reminded us that mass de- 
fect corrections are important in the degenerate plasmas 
at this level of approximation that we are interested in 
(see Blaes et al. 1992). Fortunately this is an easy fix 
as our derivation of equation ([5]) is independent of any 
assumption regarding A and Z of the ions so including 
mass defect corrections merely requires including the ap- 
propriate atomic masses of the ions0 However, we will 
show below that Coulomb corrections ultimately domi- 
nate over both thermal and mass defect corrections. 

The alert reader may worry that these mass defect 
corrections may change the results of paper I, II, and 
III. However, he or she may rest assured that these cor- 
rections are orders of magnitude smaller than the order 
unity corrections of H on C which we discussed in paper I, 
II, and III. The only possible exception may be in regard 
to our discussion of the He buffer in paper III. However 
as we will show below, Coulomb effects continue to dom- 
inate, which argues that our results on the He buffer in 
paper III remains valid. 

Figure [1] plots the He (dashed lines) and C (solid lines) 
as a function of column depth, y = J pdz, for an equi- 

5 We have used the atomic database compiled by Na- 
tional Institute of Standards and Tec hnology available at 
http: / /www. nist.gov/physlab/data/comp.cfm 



librium atmosphere with temperature, T = 10 6 K. Note 
that y = P/g in a plane parallel, constant gravity at- 
mosphere. In thin lines, we show the He and C abun- 
dance accounting for only thermal corrections to the elec- 
tric field (thin lines marked by thermal) and accounting 
for both thermal and Coulomb corrections (thin lines 
marked by Coulomb; see H3.2I for a discussion on how 
these Coulomb corrections are calculated). In thick lines, 
we also include the mass defect corrections to the masses 
of He and C. 

The analytic power law (eq. [5] for 7 = 4/3 (dot- 
ted line marked by y^ 1 ^ 2 ) matches that given by a full 
numerical calculation for He for purely thermal correc- 
tions (thin dashed line marked thermal) . Including mass 
defect corrections results in a significant change (thick 
dashed line marked thermal). The mass defect correc- 
tion drives the abundance of He away from both the 
analytic power law and numerical calculation for purely 
thermal corrections (dotted line marked by y^ 1 ^ 2 and 
thin dashed line marked thermal, respectively). In fact, 
the abundance of He increases rather than decreasing at 
large column. To understand this, we note that the mass 
defect correction for He on C is 



AF, 



defect 



m p g 



(A Ke - 2Z B 



0.003, 



(9) 



where An e is the corrected mass number of He. Thermal 
corrections, using equation ([5]), are 



Th 



m p g 



= Z V 



A c n c k B T 
'^c~ iPe 



0.002, 



(10) 



when the plasma is highly degenerate. The mass defect 
correction is fixed relative to m p g for increasing column 
because it depends only on the mass numbers of He and 
C. On the other hand, the thermal correction scale like 
ncA:BT/P e , which decreases for increasing column. Mor- 
ever, the nucleons in He are heavier than the nucleons 
in C, so the larger mass defect correction opposes the 
thermal correction, reversing the He abundance trend at 
large depth. 

Coulomb effects, however, overwhelm both the thermal 
and the mass defect corrections, which we also show in 
Figure [1] in the lines marked Coulomb. A quick order of 
magnitude estimate of the size of this Coulomb correction 
is (see the discussion in E)3. 21 especially around eqn. [33] 
for details) 



APcoi 
m p g 



3Zh c 



2 /a e 



10(7-1) E F 



0.006 



(11) 



for y — 10 10 gcm -3 and g = 2 A3 x 10 14 cms -2 and as- 
suming that nc 3> nue- Clearly Coulomb corrections 
dominate over both the thermal and mass defect correc- 
tions as the effects of these forces on the number density 
of ions is exponential. However, we have included all 
these effects in subsequent calculations in the remainder 
of this paper. 

Figured] also shows that the fraction of the He trace de- 
clines against the heavier C background, but the rate at 
which this fraction declines strongly depends on whether 
the corrections to the electric field include thermal ef- 
fects, mass defect effects, or Coulomb effects. The dot- 
dashed line is for an abundance that varies inversely with 
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Figure 1. Fractional abundance of He (dashed lines) and C (solid 
lines) as a function of y in an equilibrium atmosphere with T e = 10 6 
K and g = 2 A3 X 10 14 cms" 2 . The thin (thick) lines labeled ther- 
mal only include thermal (and mass defect) corrections and the thin 
(thick) lines labeled Coulomb contain both thermal and Coulomb 
corrections (thick lines also include mass defect corrections), al- 
though Coulomb effects dominate over both thermal effects and 
mass defect corrections. The dotted line is a y _1//2 power law 
(eq. [8]) for a relativistic electron equation of state (7 = 4/3) 
and this analytic result tracks the exact numerical calculation for 
purely thermal corrections. The dot-dashed line denotes a y~ l 
power law. The respective He columns, 3/H0 = / PUcdz, where 
PHc is the local He density, are yn c es 4 X 10 10 gcm — 2 and ~ 
7 X 10 11 gem - 2 for purely thermal corrections and for thermal and 
mass defect corrections, respectively and j/ho ~ 1-2 X 10 7 gcm — 2 
and ss 1.5 X 10 7 gcm for thermal + Coulomb corrections and 
including mass defect corrections, respectively. 

column, i.e., Hj/ntot fx For He abundances that de- 
cline shallower that this, the total He column yn e , rises 
substantially with increasing depth. For He abundances 
that decline steeper than ?/He is capped and does 
not increase substantially with increasing y. If we only 
include thermal corrections to the electric field, the He 
abundance follows the power law (eq. [8]), which is shal- 
lower than y^ 1 . Hence, j/hc continually rises with in- 
creasing y. Adding the mass defect effect to the thermal 
correction only make j/hc increase faster. However, if we 
include Coulomb physics, we find that at sufficient depth 
(y Vcut ~ 4 x 10 9 gcm~ 2 ), the abundance is steeper 
than and hence the He column does not increase 
substantially for increased y. Below y < y cu t, however, 
the He abundance is shallower than this and so the He 
column continuously increases. Hence, we find that the 
bulk of the He sits at a depth of y cn t, which we show ex- 
plicitly below. This will have important implications on 
the rate of He DNB as we show in SJU 

3.2. Coulomb Effects via a Chemical Equilibrium 
Approach 

Now we turn our attention to the effect of Coulomb 
physics on the electric field. While it would be tempting 
to follow the derivation leading to equation ([5]) , but sub- 
stitute Coulomb correction terms for the ion pressure, 
this gives erroneous results in the general case. Instead 
we have calculated these corrections starting from the 
condition of chemical equilibrium. De Blasio (2000) per- 



formed a similar calculation to study the interface be- 
tween two dissimilar ions near the crust of a neutron 
star. Our calculation below is similar in spirit, but we 
make contact between the various methods of calculat- 
ing the corrections to the electric field, i.e., starting from 
the hydrostatic balance for each ion (eq.pQ) vs. chemi- 
cal equilibrium. The condition of chemical equilibrium 
is equivalent to the equation of hydrostatic balance for 
each species, but has the advantage that the effects of 
Coulomb physics is more naturally incorporated. 

We begin by considering the free energy of an ideal 
plasma, which is 



F = F£L(N e ,V,T) + F® 



nondcj 



.(N h V,T), (12) 



where i*id,d cg and F> d nondcg are the ideal free energies 
of the degenerate electrons and nondegenerate ions, and 
N e and 2Vj are the total number of electrons and ions 
in a volume, V. For a nondegenerate ideal gas, F — 
-Nk B T(hx(V/N)^) + I), where A, = Jm^T is 

the thermal wavelength. From the free energy, we find 
the chemical potential via 



M = 



dF 



(13) 



,T,V 



Thus, the chemical potential for a nondegenerate ideal 

gas is /iW = — k B T\n (nAf), where n = N/V is the 
number density. In the presence of external fields, the 
chemical potential becomes fi = //y + <I>, where <& is the 
potential associated with the external field (Landau & 
Lifshitz 1980). For gravitational and electric fields, the 
external potential is 



$ = z{A l m p g - Z t eE), 



(14) 



where, again z — r~r rc f. Imposing chemical equilibrium 
d[ifdz = 0, we find for an isothermal, nondegenerate, 
plasma, 



kT 



d In m 



-Aiiripg + ZieE, (15) 

UT 

kT — - — - = eE, 16 

or 

which recovers the equation of hydrostatic balance for 
each species (in the isothermal limit) and demonstates 
the relationship between chemical equilibrium and hy- 
drostatic balance. 

To consider Coulomb physics, we include an additional 
term to the free energy: 



F = F$ de£ (N e , V, T)+F^ nonde£ (N i , V, T)+F ex (N e7 N l , V, 

(17) 

where F ex (N e: Ni 7 V,T) is the excess free energy due 
to Coulomb interactions. Chabrier & Potekhin (1998) 
give an excellent discussion of the Coulomb physics 
of electron-ion plasmas and we adopt their formalism 
and notation^ Following Chabrier & Potekhin (1998), 
we rewrite the excess free energy F ex = N e f ee ksT + 
Ni (fa + fie) k B T in terms of dimensionless functions, 

6 Chabrier & Potekhin (1998)'s results are suitable at this level 
of approximation, but see Potekhin & Chabrier 2000 and Potekhin, 
Chabrier, & Rogers 2009 for more accurate results. 
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/ ee (r e ), fu(T e ), and fa e (T e ), which denote the electron- 
electron, ion-ion and screening interactions respectively. 
We define T e = e 2 /a e kBT as the electron coupling 
parameter, and a e = (47rn e /3) _1 ' 3 as the mean elec- 
tron spacing. The ion coupling parameter Ti = 

Z 2 e 2 /a l k B T is related to T e via T, = Z^ /3 T e , where 

0< = (4vrn l /3)- 1 / 3 . 

Typically the ion- ion term (fa) dominates and for the 
purposes of simplifying the discussion, while capturing 
the relevant physics, we consider just this term. The 
chemical potentials then become 



dF 



(e) 



Ni,V,T 



k B T x Nj 
3 ^ N e 



u, . 



/' 



dP 



(18) 
(19) 



kBTfi, 

N e ,N j9 n,V,T 

where we have defined fa — fu(Ti) and u, = 
(r, ) /a In r 4 . Ghabrier & Potekhin (1998) (also see 
Potekhin et al 2009) have provided fitted forms for ui 
(their eq. [27]) which we use below. 

Equilibrium requires dfi/dr — ZeE — Am p g for each 
species. Ignoring temperature gradients, we find 

91nn, Uidh\n e ZieE — Aim p g 



di 



dfj,. 



d\nn, 
d In rii 



dr 

d In n„ 



k B T 



(20) 



dr 



T 

d In n c 
dr 



dr 



+ 



n e 



= -eE, 



(21) 



where hi = d\nui/d\uT. For a single ion species, n e 



and therefore d\nn 
(l2~H simplify to give 

9 Inn 




= d\nn e . So equations ([20 




1 , ^Uih 
-k B T— 
9 Z 



Am p g, 



= -eE. 



(22) 
(23) 



Including Coulomb physics introduces corrections to 
both the ion equation (|22|) and the electron equation 
(j2"3"|) . The additional correction to the electron equation 
would have been missed had we tried to derive the elec- 
tron equation from the hydrostatic equation ([2])Q The 
total pressure is P = — (dF/dV)\ T . Therefore, the non- 
ideal, Coulomb part is Pcoul = nkTui/3, which agrees 
with the sum of equations (|22|) and (|23| . 

This is further complicated in the general case of mul- 
tiple species where the relation d In rii — d In n e does not 
hold. For instance, consider a trace ion species against a 
background Zim n e , equation (|20j) becomes 



9P„ 



(id) 



dr 



A 



im p g 



Z,eE 



1 51nn e 
-k B Tui— - — 
3 or 



where P. 



(id) 



(24) 

riik B T is the ideal gas equation of state. 
Here a simple substitution of a corrected ions pressure, 

7 The ion equation J22H is also different from what we would 
have derived had we started from the ion hydrostatic equation (TTJ . 



Pi = riik B T(l + Ui/3), would have given erroneous re- 
sults. Distinct terms must be included for both the ion 
and electron momentum equations. The Coulomb inter- 
action must not be regarded as a pressure for each ion 
species, but rather as a force that is sensitive to the elec- 
tron gradient which is set by the background. 

This discussion of these two simple cases (single ion 
species and trace) highlights the problems with calculat- 
ing the electric field starting from the equations of hydro- 
static balance for ions and electrons (eq. [T] and [5] ) when 
Coulomb physics plays a role. Fortunately, our approach 
from the viewpoint of chemical equilibrium avoids these 
problems. 

We now solve the general case by first plugging the ion 
equation into the Coulomb term of the electron equa- 
tion. Collecting terms, the result can be written in the 
intuitive form 



d In '/ 



ZZeE 



lm p g 



E* 



where 



e: 



a: 



nsUiZi 



d\nn, 

1 X ^ 

~3^ a 

i 

_1 rijUtAj 



(25) 



- 1 ,(26) 



(27) 
(28) 



This result can now be used to simplify the ion equation: 
d In rii _ Z*eE — A*m p g 



where 



dr 



Z* = Zi - 



A* = Ai 



k R T 



k B T UjZ* 
Et 3~ 
k B T u t A* e 
~E*. 3 



(29) 

(30) 
(31) 



In this form, it is as if the charge and mass of each species 
(ions and electrons) are shifted due to Coulomb interac- 
tions. Using charge neutrality, we solve for the electric 
field and find 



eE Y Ji ^i/n e )Z i A* - (k B T/E* e )A* e 



m P g Y.Mi/n^ZiZf - (k B T/E*)Z* e 



(32) 



One can easily show this equation has the cor- 
rect limits for noninteracting particles, for degener- 
ate/nondegenerate and nonrelativistic/relativistic elec- 
trons. 

As an illustrative exercise, we use equation (|32[) to find 
the electric field for a single ion species: 



eE 



A[(j- l)E F + k B Tuh/9Z] 



m p g Z [(7 - l)E F + k B Tuh/9Z] + k B T (1 + u/3) 



A 
Z 



where Pjf oc n e Ep oc nj is the degenerate electron 
pressure. In analogy with the thermal correction to the 



3 Z 2 / 3 e 2 /a ( 
10(7 - 1) E F 



(33) 
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electric field, the correction to the electric field due to 
Coulomb physics scales as n,iEc/n e Ep, where Eq is the 
Coulomb energy of the ions. 

Equation (|52|) is the central result of this section 
and we now briefly discuss its essential physics. When 
Ai/Z\ ^ A^jZi-, i.e., dissimilar p e , the ion scale height 
is determined by ksT /m p g which is much smaller than 
the electron (pressure) scale height l e — Ep/m p g in 
the strongly degenerate part of the atmosphere. When 
Ai/Z\ sa A2/Z2, i.e., similar p e , then the zeroth or- 
der field produces no differentiating force (modulo ther- 
mal/mass defect corrections). Moreover, the thermal 
correction to the electric field produces an ion scale 
height w l e as shown earlier. Mass defections correc- 
tions produce ion scale heights that are similarly large. 
While this results in a decreasing fraction for the lighter 
ion going deeper into the envelope, the column of this 
lighter ion increases as discussed earlier (and increase re- 
ally fast if we include mass defect corrections). Including 
Coulomb physics causes the ion scale height to become 
l e /Ti, much smaller than l e for large T, but is still larger 
than the ion scale height for dissimilar A/Z by Ec/Ep. 
Both concentration and number density for the lighter 
ion decrease with increasing depth. 

We now utilize the discussion above to estimate j/ C ut- 
Plugging equation (|3"3")l into hydrostatic balance (JT|) for 
a trace, 2, and a background, 1, we find for relativistic 
degenerate electrons (7 = 4/3) that 



<91nrai 

dr E F 
d In 712 rn p g 



m p9 r 7 5/3 



dr 



E F 



5/3 



(34) 
(35) 



Subtracting (|3"4")l from ([33]) and noting that d In P/dr — 
dlny/dr = -pg/P = -4Am p g/ZE F B we find 



d In nijnx 
dhiy 



1 Z 



Z. 



5/3 



zi 



5/3 



(36) 



Plugging in values for a He/C envelope appropriate for 
Figure [T] with base temperature Tb = 1 x 10 s K (see £Q] 
for a discussion of the base temperature), we find that 
dln(n2/ni)/dln y = — 1 at p = p cut ~ 5 x 10 6 gcm~ 3 
for relativistic degenerate electrons, which gives y cut w 
2xl0 9 gcm- 2 . 

The results of this section show that stratification oc- 
curs even for elements with the same p e . This has impor- 
tant implications in high gravity environments like white 
dwarfs and NSs beyond the problem of He DNB. 

4. DIFFUSIVE NUCLEAR BURNING OF HELIUM 

Combining the results of the thermal profile of the 
envelope, rates of He burning, and rates of diffusion, we 
now calculate the rate of nuclear burning of He in NS 
envelopes. For the thermal profile, we solve the heat dif- 
fusion equation for a constant flux envelope in radiative 
equilibrium: 

OT 3np ^ (3?) 



dr 



16T 3" e> 



8 The factor of 4 comes from noting that P fs P e , i.e., relativistic 
degenerate electrons dominate the pressure, and P e = n e Ep/A. 



where n is the opacity. Starting from the photosphere 
(with optical depth, t = 2/3), we integrate the plane- 
parallel envelope inward assuming constant flux to a 
depth of y = yt = 10 14 gcm -2 , which we call the base 
of the envelope. As in paper II, we use the tabulated 
conductivities of Potekhin et al. (1999) and the radia- 
tive opacities and the equation of state of Potekhin & 
Yakovlev (2001), which is applicable to all magnetic field 
strengths (see paper II for more details). 

Potekhin ct al. (1997) performed a similar calculation 
down to a density of p = pb = 10 10 gcm~ 3 , which they 
referred to as the base of the envelope. These two defini- 
tions are roughly equivalent to each other, namely that 
p{Ub) ~ 10 10 gcm -3 . We define the temperature at g/& 
as the base temperature, Tf,. For a NS in thermal equi- 
librium, the base temperature and the core temperature, 
T c are equal to each other as the temperature profile is 
nearly isothermal due to the large conductivities. For 
the most part, Tb and T c can be use interchangable, ex- 
cept for the very early cooling history of the NS (age 
< 100 yrs), when the NS is not in thermal equilibrium 
and Tb ^ T c . This point will be important in our discus- 
sion of Figure [5j 

For the diffusion coefficient of the trace He, we use the 
results of paper II (also see Brown, Bildsten & Chang 
2002), where we found: 



V « 10 



-3 



A^Tj 3 

vl.3 70.3„0.6 
Z l Z 2 Pb 



cm 2 sec , 



(38) 



where T 6 = T/10 6 K and p 5 = p/10 5 gcnr 3 . As dis- 
cussed in paper II, equation (|38j) is reasonably accurate 
in the liquid regime of the background C, 1 < T < Tm, 
where = 175 defines the melting point of the crys- 
talline phase (Potekhin & Chabrier 2000). For T > T M , 
the material is a crystalline and we assume that there 
is no diffusion, though we never reach this point in our 
parameter regime. 

We now discuss the nuclear processes that consume 
He. For our nuclear reaction rates, we have utilized both 
the NACRlfiAngulo et al. 1999) and REACLIl£3 com- 
pilation of nuclear reaction rates and experimental val- 
ues (see nucastro.org). We have not included electron 
screening in our calculations, though the effect of screen- 
ing would be to increase the rate of burning. For triplc- 
a reactions, we use the fit in Nomoto, Thielemann, & 
Miyaji (1985), but have found this rate to be typically 
very small for the parameters of interest because of the 
reduced number density of He in the burning region (note 
the reduced number fraction of He in Figure [1]). Hence, 
we no longer discuss triple-a reactions. 

The local rate of He burning is 



Z/Hc 



Z/Hc 



THccol 



dz- 



4n He TOp 



) ^Q — Cap; 



(39) 



where z is the depth and n Q _ cap is the density of a He- 
capturing substrate. 

In the top panels of Figures [2] and [3l we plot the rate 
of He burning (solid lines) as a function of column, i.e., 
dyne/dy for effective temperatures of T e = 1.25 x 10 6 K 
and 2 x 10 6 K with photospheric He abundances of 50%. 



9 http:/ /pntpm. ulb.ac.be/Nacre/nacre. htm 

10 http:/ /nucastro.org/reaclib.html 
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We also plot the yn e (short-dashed lines) as a function of 
column. Plotted in the bottom panels of these two fig- 
ures is the background temperature profile (dotted lines) 
and the composition profile (long-dashed lines) of the 
envelope. Figure [2] and [3] represent He DNB in the nu- 
clear limited regime and the diffusion limited regime re- 
spectively (see 35] and paper II for a detailed discussion 
of the nuclear limited and diffusion limited regimes of 
DNB). In both cases, the peak of the burning is set 
by a combination of a rising temperature profile and 
a rapidly falling He concentration. The cutoff in the 
nuclear-limited regime (Fig. [5]) results from the equilib- 
rium profile of He embedded in a C background. Hence, 
the He profile remains in diffusive equilibrium as it is 
slowly consumed. 

On the other hand, the sharp dropoff in the He pro- 
file in the diffusion limited regime (FigJ3]) is not a result 
of the equilibrium structure, but rather due to nuclear 
burning forcing the He profile out of diffusive equilib- 
rium. As in the case of diffusion limited H-DNB dis- 
cussed in paper II, the rate at which H is consumed is 
ultimately limited by the rate at which H diffuses down 
to the burning layer. Below this layer H does not exist as 
the burning layer captures all downwardly diffusing pro- 
tons. The rate of nuclear burning is not limited by the 
rate of nuclear reactions, but rather by the rate of diffu- 
sion down to this boundary. We see this same behavior 
for He at sufficiently high temperatures. 

To calculate the case of diffusion limited He DNB, we 
have generalized the argument of paper II. The details 
can be found in paper II, but for completeness, we write 
the additional equation to be solved below: 



d 2 / Hc 
dz 2 



dfuc 9 In n Hc 
dz dz 



/hc 



(40) 



where riHc.o is the He number density in the absence 
of nuclear reactions, /n e = ^He/«He,o is the correction 
factor in the presence of nuclear reactions, and njjc is 
the Hc number density. In the regime where equation 
(|4"D|) becomes important, the second term on the LHS 
becomes small compared to all the other terms. In other 
words, the scale height of /ne is small compared to the 
equilibrium scale height, h cq — (dhin^e.o/dz)^ 1 . As 



we found in paper II, equation (|40|) becomes a diffusion 
equation in the presence of a nuclear driven source. 

We note that the small difference between the effective 
temperature between Figure [2] and [3] (1.25 x 10 6 K vs 
2 x 10 6 K) yields a seven order of magnitude difference in 
burning rate and He column lifetime. This results from 
the strong temperature sensitivity of a capture onto C. 

We now discuss how the lifetime of a He layer, th c = 
J/He/^He) scales with its column, yn c . In Table [TJ we 
show the photospheric He abundance (nne/^tot at the 
photosphere), He column yn c (for He sitting on top of C), 
and the He lifetime, t Hc , for a NS with T e = 2 x 10 6 K. We 
note two remarkable things about this table. First, the 
He lifetime is independent of the total column. This was 
not the case for H DNB (papers I, II, III). Second, even 
though the photospheric abundance of He can be low, the 
total column of He is high. Namely, for a photospheric He 
abundance of 10%, the total He column is ~ 2500 gem" 2 
much greater than the photospheric column of y p h » 
1 g cm" 2 . 




5 10 
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Figure 2. Profile of He burning (solid line) and j/jjc (short-dashed 
line) as a function of column with T e = 1.25 X 10 6 K (upper panel). 
Shown in the lower panel are temperature (dotted line) and He 
abundance (long-dashed line). The base temperature is 7.8 X 10 7 
K, photospheric He abundance is 50%, and the total He column 
depth is jHc ~ 2100 g cm -2 . At this base temperature, the rate of 
He burning is limited by the rate of alpha capture onto carbon, i.e., 
the nuclear limited regime of DNB. The He tail which penetrates 
into the C remains in diffusive equilibrium as the rate of alpha 
capture is slow. 
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Figure 3. Same as Figure [2] but with T e = 2 X 10 6 K. The base 
temperature is = 1.7 X 10 8 K, photospheric He abundance is 
50%, and the total He column depth is jh c ~ 4600 gem -2 . Unlike 
Figure [2] the rate of He burning is limited not by the rate of alpha 
capture onto carbon, but on the rate of He diffusion down to the 
burning layer, i.e., the diffusion limited regime of DNB. The He tail 
which penetrates into the C is not in diffusive equilibrium below 
the burning layer. 

We now demonstrate that the independence of Tn e and 
2/He and yn c 3> j/ p h are due to the same physics. The rea- 
son why ?/He 3> y P h is because the bulk of the He is not 
near the photosphere, but rather deep in the envelope. 
We show this in the top panels of Figures [5] and [3] where 
we plot ?/He (short-dashed lines) as a function of y. Note 
that the bulk of the He does not reside near the sur- 
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Table 1 

Photospheric He abundance, He column, and corresponding 
lifetime for a He/C envelope on a NS with T e = 2 X 10 6 K. 



He photospheric abundance 


lgio2/He [gem 2 ] 


IgloT&e [y rs l 


0.1 


2.888 


1.947 


0.3 


3.397 


1.947 


0.5 


3.664 


1.947 


0.9 


4.144 


1.947 


0.999 


4.851 


1.947 


0.999999 


5.848 


1.949 



face, but increases substantially until is asymptotes at a 
column of y « y cut . 

We now explain why the bulk of He resides near y cu t- 
In paper I, we showed that the total lifetime of H on a 
proton capturing substrate is dominated by the time it 
takes to remove a photospheric H column. Namely, we 
found that the H lifetime, tjj, scales with the H column 
as th oc y^~ S , where 6 = A-x{Z\ + l)/A\ — Z-x — 1, where 
1 and 2 denote the background and trace ion species re- 
spectively, i.e., 1 + 5 = —5/12 < for a H on C envelope 
(paper I). As we discussed in §5.2 of paper I, this scaling 
arises from the power-law falloff of H abundance in a non- 
degenerate atmosphere, which is fT-jr/i-tot oc y s = y~ 17 / 12 
(see eq. [31] of paper I). For He on 0, the He abundance 
follows nue/fhot oc y s — y~ 2 ^ 3 . The H abundance on a 
C substrate falls off steeper than y" 1 , whereas the He 
abundance on a C substrate falls off shallower than y -1 . 
Thus, the H column is capped at the H/C boundary, the 
column where the number density of H and C are equal, 
but the He column continues to increase substantially 
below the He/C boundary. Indeed, one can show that 
the He column always increases below the He/substrate 
boundary for any substrate with A/Z w 2. This increas- 
ing He column carries through to the Coulomb case as 
discussed earlier in t j3.2l Indeed, the He column contin- 
ues to rise until it is capped at a column of y cu t as we 
argued in i j3.ll 

We now note that the burning layer is close to the re- 
gion where He begins to get cut off, j/bum ~ 2/cut (see 
for instance, Figures [5] and ■ If we say that they are 
the same, we are left with the result that the He abun- 
dance in the burning layer and the total He column are 
linearly related to each over. Since the rate of He DNB 
(via captures onto C) is linearly related to the He abun- 
dance (and hence yn e oc yHe), we find that the lifetime 
of the He layer, th = ync/yne is independent of yne as 
demonstrated in Table [TJ explaining the independence of 
r Hc and y Hc - 

In Figure^] we show the He lifetime, th , for a range of 
a-capturing material as a function of base temperature, 
Tj, (left), effective temperature, T e , and 5 = (solid 
lines). We also plot (dotted lines) the same case but 
for a radial B-field of B = 10 12 G. Magnetic fields up 
to pulsar strengths do not make a significant difference. 
For sufficently large base temperatures, the lifetime to 
deplete the atmosphere of He is a few weeks. 

To compare the lifetime of a He layer with the typical 
cooling history of a NS, we plot in Figure [S] a few repre- 
sentative cooling tracks (long-dashed lines) and overlay 
the THe-7b relation shown in Figure[4j Here, we have plot- 
ted Tf,, which during the very early history (age < 100 



yrs) of the NS ^ T c . Only after the cooling wave sweeps 
through the star, which can be seen by the sudden drop 
in TJ, at ~ 100 yrs, is the NS in thermal equilibrium. As 
we have done in paper II, we include a standard cooling 
track which presumes modified URCA cooling on a 1.3 
Mq NS, another with core proton superfluidity (Potekhin 
et al. 2003), and one with triplet-state neutron superflu- 
idity in the core with a maximum critical temperature of 
8 x 10 8 K (calculated by D. G. Yakovlev previously for 
paper II). 

The remarkable aspect of this plot is that independent 
of the cooling model assumed (where we have taken a few 
representative models that span the range from standard 
cooling to fast cooling), the lifetime of a He atmosphere 
is short compared to the cooling age (during its early 
history). For instance, He on a C envelope will be de- 
pleted even in the fast cooling case, if it was deposited up 
to 1000 yrs after formation. For slowly cooling models, 
this can extend up to almost 1 Myr. For a much heavier 
material, i.e., Ar, the lifetime is still short ~ 100 yrs for 
both fast and slow cooling. Thus, primordial He should 
be depleted on NS atmospheres during it early cooling 
history. 

As we close this section, we note that the lack of ob- 
served light elements on NS photosphere may not indi- 
cate that there are no light elements present. As we 
discussed above, in diffusive equilibrium, the bulk of the 
He resides near y cu t- What this implies is that while the 
photospheric abundance of He could be low, the total 
He column could be large (see Table [IJ. This result has 
implications for the subsequent accretion or production 
of light elements on NS surfaces. For instance, suppose 
a sudden accretion or spallation event produces a total 
He column of yn e ~ 10 3 g cm -2 on the NS surface, which 
previously consisted of C. The NS surface would initially 
appear to be composed of light elements. Over time, 
as the He diffuses downward to y cu t, the photospheric 
abundance of He drops and the NS surface appears more 
C-like. 

The timescale for this mixing is the diffusion time down 
to y cu t for He on C: 

h 2 

"p.cut » 1.3^1.3 Tr _„ ( A -,\ 

Tditf ,cut ~ ~j ~ 2 p 7 l s yis, (41) 

where /i p , C ut is the pressure scale height at y cu t, Pi — 
p/10 7 gcm- 3 , and T s = T/10 8 K. This timescale is en- 
ticingly similar to the timescale of the observed spectral 
variations in RX J0720.4-3125 (de Vries et al. 2004; Vink 
et al. 2004; Haberl et al. 2006; Hohle et al. 2009). 
Haberl et al. (2006) proposed that NS precession may 
be responsible for RX J0720.4-3125's observed spectral 
hardening and subsequent softening. However, we specu- 
late that slow evolution of the NS surface over a period of 
years may be the culprit. Hydrogen on RX J0720. 4-3125 
is consumed by DNB, whereas He would diffusively mix 
down to y cut . A more detailed study of this process of 
diffusive mixing for He is needed to make a more quanti- 
tative comparions, but such a calculation is well beyond 
the scope of this paper. 

5. CONCLUSIONS 

We have now extended our prior work on H DNB to 
He. We now include thermal, mass defect, and Coulomb 
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Figure 4. Lifetime for He sitting on top of C, Si, or Ar envelope as a function of base temperature, T b (left plot), and as a function of 
effective temperature, T e (right plot) for B = (solid lines) and B = 10 12 G (dotted lines). For sufficiently high T b (T e ), i.e., T b > 1.6 X 10 8 
K (T e > 2 X 10 6 K), the lifetime of He against capture on C is very short (< 10 2 yrs). As these Tj/s are reached during the early cooling 
history of the NS, primordial He should be burned off. At low temperatures, there is a steep power-law dependence, as would be expected 
for DNB in the nuclear-limited regime. At higher temperature, the dependence on T b becomes shallower, indicative of DNB in the diffusion 
limited regime. 

the absence of external sources. In addition, the effects 
of magnetic field, as discussed in (also see papers II 
and III for H DNB), do little to change this basic result. 
We require the presence of H and He capturing material 
in NS envelopes, which the initial nuclear evolution of a 
cooling NS produces (Hoffman & Heyl 2009), i.e., proton 
and He capturing elements as light as Si. Finally, we 
comment on the timescale for He on C to reach diffusive 
equilibrium and suggest that such an evolution may ex- 
plain the observed long-term spectral variation seen in 
RX J0720.4-3125. 

Since H and He are consumed during the NS's early his- 
tory, we expect that the NS surfaces will be dominated 
by mid-Z elements such as C, N, and O. The evidence 
for mid-Z elemental compositions on NSs is sparse. Ob- 
servations of 1E1207.4-5209 by Sanwal et al (2002) and 
Mereghetti et al. (2002) found X-ray spectral lines at 
0.7 keV and 1.4 keV with additional spectral features at 
2.1 and 2.8 keV, which have been disputed (Mori et al. 
2005). Subsequent modeling by Hailey & Mori (2002), 
Mori & Hailey (2006), and Mori k Ho (2007) suggested 
that its atmosphere consists of mid-Z elements like O or 
Ne with a magnetic field strength of6x 10 11 G. Given 
the age of 1E1207-52 of » 7 kyrs (Roger et al. 1988), 
we found that H is easily consumed (paper II). Using 
the results of this paper, we suggest that He could have 
been consumed during the early cooling history of the 
NS exposing the underlying O or Ne. Subsequent accre- 
tion was O-rich and not subject to spallation or that the 
strong pulsar wind prevents accretion. We also note that 
removal via pulsar winds remains a strong possibility. 

However, more recent work by Gotthelf & Halpern 
(2007) challenges the intepretation that these spectral 
lines on 1E1207-52 are due to mid-Z materials. Rather 
the spindown rate which Gotthelf & Halpern (2007) 
found suggest a B-field strength (B < 3 x 10 11 G), con- 
sistent with electron cyclotron lines for the higher energy 
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Figure 5. Base temperature as a function of He lifetime/age. 
Overlayed are cooling tracks (long-dashed lines) for various NS 
models (Potekhin et al. 2003) including standard cooling (no su- 
perfluidity, modified Urea cooling represented by 1.3iWg noSF), 
proton superfluidity (pSF) and triplet-state core neutron superflu- 
idity (nSF). Note that for sufficiently young NS or sufficiently slow 
cooling, thc < age, which indicates that any He on the NS will be 
consumed by DNB. 

corrections to the electric field to calculate the structure 
of He on C (or any other « 2 material) in diffusive 
equilibrium. We show that the bulk of He in diffusive 
equilibrium sits at a layer y cut and show that as result 
of this the lifetime of a He on a NS is independent of 
its column. We demonstrate that young NSs undergo 
a phase where He is consumed via DNB (see especially 
Figure [S]). Combined with the results of paper I, II, and 
III (see Figure 8 of paper II), we conclude that all pri- 
mordial H and He are depleted on the surfaces of NS in 
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(1.4 keV) line. 

Currently the best evidence for mid-Z atmospheres 
comes from Ho & Heinke (2009). They recently argued 
that the surface of the NS in the center of Cassiopeia A 
consist of a C atmosphere. Their modeling of archival 
observations of this compact X-ray source showed that 
a T c = 1.8 x 10 6 K carbon atmosphere NS with a low 
magnetic field provides both a good fit to the observed 
spectrum and is consistent with the theoretical expecta- 
tion of the radii of NSs, i.e., R = 12 — 14 km. Additional 
observations of other young NS may provide further ev- 
idence of mid-Z elements in their atmospheres. 

We should mention that even though DNB should de- 
plete the atmospheres of young NSs of any primordial 
H or He, subsequent accretion with or without spalla- 
tion (Bildsten, Salpeter, & Wasserman 1992) could still 
lead to a H/He atmosphere in spite of DNB. As the mass 
of the photosphere is low and the timescales still fairly 
long, even a small amount of accretion will pollute the 
photosphere. A quick estimate of the amount of accre- 
tion needed to overwhelm DNB is 

M > 10- 30 ( Vph A ( -^-) Mq yr-\ (42) 
Vlgcm- 2 J VlOOyrs J ° y v ' 

A pulsar wind from the NS may suppress this accretion, 
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but the amount of suppression has to be fairly strong as 
indicated by equation (|42|) . 

Armed with the results of this work and prior efforts, 
we expect that primordial H and He are depleted, and so 
any observed H or He on the surfaces of these NS must be 
due to subsequent accretion (with or without spallation). 
If this subsequent accretion can be stopped or prevented 
via a pulsar wind, for instance, the underlying mid-Z 
material would be exposed. 
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